Sub-diffraction light propagation in fibers with anisotropic dielectric cores 
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We present a detailed study of light propagation in waveguides with anisotropic metamaterial 
cores. We demonstrate that in contrast to conventional optical fibers, our structures support free- 
space-like propagating modes even when the waveguide radius is much smaller than the wavelength. 
We develop analytical formalism to describe mode structure and propagation in strongly anisotropic 
systems and study the effects related to waveguide boundaries and material composition. 
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I. INTRODUCTION 



Guided propagation of optical signals in isotropic dielectric fibers is possible only when the light-transmitting region 
' of a fiber is at least of the order of the free-space wavelength^. This fact strongly limits resolution of modern optical 
microscopy and spectroscopy, prevents the construction of compact all-optical processing units and further develop- 
ment in other areas of photonics 2 *^^^. Although it is possible to propagate GHz radiation in deep subwavelength 
i - 1 \ areas in coaxial cables or anisotropic magnetic systems&^ilSiil, the direct scale-down of these techniques to optical or 
q • IR domains is problematic^. 

• i-H ' Until recently, all designs involving optical light transport in subwavelength areas relied on the excitation of a surface 
O ., wave - a special kind of electromagnetic wave propagating at the boundary between the materials with positive and 

negative dielectric constantsi^i^iSii&. The spatial structure of surface waves, however, is fundamentally different 
, from the one of "volume" fiber modes or free-space radiation. While it is possible to couple the radiation between 
the volume and surface modes, such a coupling is typically associated with substantial scattering losses and involves 

* c/2 \ materials and devices of substantial (in comparison with optical wavelength) siz o 1 ! 1 ? , 

A new approach to guide volume modes in subwavelength areas has been recently introduced in Rcf. 1 ^. It has been 
demonstrated that the waveguide with strongly anisotropic dielectric core supports propagating modes even when its 
radius is much smaller than the wavelength. The phase velocity of these propagating modes can be either positive or 
negative depending on the waveguide core material. It has been also shown that anisotropic core waveguides can be 
t-H , tapered and used to effectively transfer free-space radiation to and from nanoscale. These tapered systems, photonic 
^ ■ funnels, are in a sense variable-index (n 0), volume mode-analog of the adiabatic surface mode compressor, proposed 
; in Refi^. 

In this manuscript we describe the perspectives of light propagation in anisotropic waveguides. In Section 2 we 
| describe the physics behind the volume mode propagation in subwavelength waveguides with perfectly conducting 
metallic walls. Effects of the core micro-structure on the light propagation are considered in Section 3. Section 4 
\Q ' is devoted to the effect of waveguide walls on mode structure, confinement, and propagation constants. Section 5 
f^) , concludes the paper. 



II. LIGHT TRANSMISSION IN WAVEGUIDES WITH PERFECTLY CONDUCTING WALLS 



C/3 1 

■ The geometry of the light transmission considered in this work is schematically shown in Fig. ^ The direction 
of light propagation is assumed to coincide with z axis of cylindrical coordinate system. The optical axis of the 
anisotropic core material is assumed to be parallel to z direction; en and ej_ denote the values of the effective dielectric 
permittivity parallel and perpendicular to optical axis respectively. 

Following the formalism, introduced iniSiiS, the waves propagating in a waveguide with perfectly conducting walls 
can be represented as a series of (TE) and (TM) modes. Each mode has its own spatial structure across the waveguide, 
determined by a structural parameter k as described below. The propagating component of the mode's wavevector 
can be related to its frequency through the free-space-like dispersion relation: 

w 2 

k 2 z =tv—, (1) 
<r 

where to is the angular frequency of light, c is the speed of light in the vacuum, and the polarization-dependent 
propagation constants 2 ^ 
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FIG. 1: Schematic geometry of a waveguide with strongly anisotropic dielectric core (a). Nano-layered metal-dielectric meta- 
material realization of a strongly anisotropic core is shown in (b). Optical axis and the direction of light propagation are parallel 
to z axis. 



Vtm = f n ■ 

For TM waves the structural parameter k and the mode structure across the waveguide can be related to the z 
component of E field, which is in turn is determined from the differential equation 

A 2 E Z + k 2 E z = 0, (3) 

with A2 being the 2D Laplacian operator, and an additional condition that the electric field E z satisfies the boundary 
conditions along the core-cladding interface 1 . Although Eq. determines mode structure in a waveguide with 
arbitrary cross-section, here we consider a cylindrical waveguide to illustrate our approach. In this geometry the 
structure of TM modes is described by E z (r, <fi, z) oc J m (K^ M r)e im ^ +jfciZ . Similar considerations relate TE waves to 
H z (r,(j),z) oc J m (K^ £ r)e lm * +lfczZ . The structural parameters are: 

o _ X(TM\TE) , s 

K (TM\TE) — ^> W 

where R is a waveguide radius and X is given by J m (^TAf) = for TM waves J' m {XTE) = for TE waves respectively. 

The dispersion equation Eq.Q is fundamentally similar to the dispersion of a plane wave in isotropic material. 
Indeed, the combination n = k z c/uj = zty/eU, which plays a role of the effective index of refraction, is combined from 
two (mode-dependent) scalar quantities, e and v. The mode propagation requires both propagation constants to be 
of the same sign. 

While the parameter e depends solely on the dielectric properties of the core material, the propagation parameter v 
can be controlled (through k), by changing the waveguide (or mode) geometry. Since k is inversely proportional to the 
waveguide size [see Eq.(0J)i, there exists a cut-off radius R CI ~ A/2, corresponding to v = for every free-space-like 
mode propagating in the waveguide with isotropic dielectric core. The modes propagate in structures with R > i? cr , 
and are reflected from thinner systems. This appearance of the cut-off radius in all dielectric waveguides can be 
considered as a manifestation of a diffraction limit - it is impossible to localize a combination of plane waves to a 
region much smaller than the wavelength inside the material A = Ao/n, with Ao being free space wavelength. 

Material anisotropy in a sense opens another direction in controlling the light propagation in waveguides. Indeed, 
anisotropy of dielectric constant makes the TM-mode parameters e and v completely independent of each other (TE 
waves are not affected by material anisotropyi&IS) . Extremely anisotropic optical materials may be used to achieve 
the volume mode propagation in deep subwavelength waveguides. 

As it directly follows from Eas. (|ll2|l . when (e± > 0, en < 0), the parameter v is positive regardless the size of 
the system. Thus, the cut-off radius does not exist at all: the decrease of the waveguide radius is accompanied by 
the decrease of the internal wavelength of the mode 2it/k z oc R, in a sense "postponing" the diffraction limit in the 
system. 

The case of opposite anisotropy (e± < 0, ey > 0), is of a special interest. The mode propagation is now possible 
only when v < 0, which in turn requires the waveguide cross-section to be extremely small. Furthermore, causality 
arguments now require the phase velocity of a propagating mode to be negative (Fig^ 181 ^. In a sense, such a 
waveguide is a complete antipode of a conventional dielectric fiber, in terms of phase propagation, as well as in terms 
of cut-off radius. 

FigElshows the dependence of refractive index on the waveguide radius. It is clearly seen that the "subdiffractional" 
light propagation indeed follows from the reduction of the internal wavelength [increase of effective n] in thinner 
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FIG. 2: Dashed lines: dispersion relation of TMoi mode in a waveguide with Si core; R = 80nm. Solid lines: dispersion relation 
of the TMoi mode in the waveguide with anisotropic core [material dispersion is neglected]; k = u/c. Material parameters are 
those of Si-Ag composite described in the text for Ao = 1.2pm (a) and Ao = 500nm (b). The former case corresponds to a 
negative refraction (ex < 0, ey > 0), while the latter describes positive refraction case (ej_ > 0, ey < 0). Positive group velocity 
is assumed. See text for details 




FIG. 3: Effect of material inhomogeneity on the mode propagation in the waveguides with homogeneous anisotropic cores 
(solid lines) and in waveguides with nanolayered composite cores (dashed lines). Material parameters correspond to TMn 
mode propagation in Si-Ag system, Ao = 1.2/xm (a) and TMoi mode propagation in Si-Ag system, Ao = 500nm (b). Horizontal 
and vertical lines correspond to k z (a + b) — 1 and to htm {a + b) = 1 respectively. 



waveguides. The figure also compares the effective index of refraction in homogeneous anisotropic material to the one 
in anisotropic meta-material composed from a number of thin metal-dielectric layers as described below. 

III. MATERIAL INHOMOGENEITY AND SUBWAVELENGTH LIGHT PROPAGATION 

Confinement of the propagating waves to deep subwavelength areas requires extremely strong anisotropy of the 
dielectric core. While a large number of anisotropic materials exist in nature, only a few of them exhibit the required 
strong anisotropjiSS. In this section we consider the perspectives of using nanostructured composites, known as meta- 
materials, as strongly anisotropic dielectric cores. Some examples of these systems include nanolayer and nanowire 
struct ur o 21 i 22 . Here we consider the case when strong anisotropy is achieved in an array of metal-dielectric nanolayers, 
shown in Fig^jD. 

To fabricate these structures, one may use standard CVD, MOCVD, or e-beam writing techniques to build a 
nanolayer composite with the total height equal to the length of the anisotropic fiber, followed by ion-beam or 
chemical etching to "shape" the cylindrical or conical waveguide. E-beam writing or self-assembling techniques could 
be also used to directly fabricate the fiber with nanolayer core. 

When the characteristic thickness of each individual layer becomes much smaller than the wavelength, the properties 
of the metal-dielectric composite are well-described by effective-medium theory (EMT) 21 : 

e ± = (5) 

ad + a m 

(ad + a m )e d e 

e \\ = ; ' 

ad^m + a m e d 

where a d , a m and e d > 0, e m < are thicknesses and permittivities of dielectric and metal layers respectively. 

As an example, in this work we assume that the core material is composed from 15 — nm-thick layers of Ag (modeled 
using Drude approach) and Si (esi = 12). As it directly follows from Eq.(|SJ), this system works as right-handed 
(ey = -59.48 + 2.78i, e x = 1.72 + 0.06?) for A = 500 nm, and as left-handed (ey = 28.72 + 0.12^, e ± = -30.51 + 0.77i) 
for Ao = 1.2 fj,m. The mode behavior in these systems is illustrated in Fig.EI Two approaches are used to calculate the 
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propagation constant of the mode in each system: EMT [Eq.JSJ], and analytic solution of light transmission though 
a ID periodic layer arrajii&Sii. It is seen that the predictions of both techniques are almost identical for the thicker 
waveguides, but strongly disagree for thinner systems. In fact, inhomogeneous microstructure of the waveguide core 
introduces the cut-off radius in anisotropy-based systems. While the appearance of such a critical radius may seem 
similar to the one in "conventional" dielectric fibers, the two have fundamentally different origins. In homogeneous 
systems the wave propagation becomes impossible when one tries to confine the propagating wave to the spatial area 
smaller than the wavelength. In metamaterial-based structures, on the other hand, the wavelength "self-adjusts" to 
the waveguide radius. The mode cut-off in this case has its origin in the break-up of the EMT when the internal field 
variation scale [A = min(27r/fc 2 , 2w/k)] becomes comparable to an inhomogeneity scale [A ~ (ad + a m )^. 

We note that while thinner layers may in principle allow unrestricted confinement of free-space light, in reality 
such a confinement will be limited by finite-size corrections to the material properties of the layers (spatial dispersion, 
Landau damping2&S£) . For metals, the minimum layer thickness can be estimated using 

a ™ ~ Ao T Too' (6) 

with Vf being Fermi velocitji2&. 

The range of wave propagation in metamaterial-core waveguides is affected not only by spatial dispersion, but also 
by implicit material absorption. For the systems that can be successfully treated with EMT the field attenuation given 
by the imaginary part of the propagation constant k z depends on waveguide geometry along with material properties 
[see Eq.Q)]. For TAfoi modes in Ag-based systems with R = 80nm considered here the attenuation is dominated 
by absorption inside Ag. We estimate that the intensity of the light attenuates in e times on the length of 1.5/xm. 
This attenuation, although it is acceptable for the short-range signal transfer, may be further compensated or even 
reversed, by implementing the gain into the "dielectric" component of a metamateriati^. 



IV. EFFECT OF WAVEGUIDE WALLS 



In this section we consider the mode dynamics in waveguides with dielectric or metallic walls. Similar to the case of 
perfectly conducting walls described above, the light propagation in fibers with any isotropic cladding can be related 
to the propagating waves with TE and TM polarizations. In this approach, the field is represented as a series of 
TM and TE waves with same propagating constant k z and frequency u>, and the boundary conditions are used to 
find the effective refractive index n — k z c/uj. Note that TE and TM components of the mode have similar but not 
identical spatial structure inside the anisotropic core. Explicitly, this structure is given by J m (K(T e\t M)r) exp(im^) 
with k^ e = e±uj 2 /c 2 — k 2 , and k^ m = t\\(u) 2 jc 2 — k 2 /e±). The mode structure in the cladding material is described 
by K m (K c \r) enp(im(/)) with k 2 x = k 2 — e c \ui 2 /c 2 , and e c i being permittivity of the cladding. The boundary condition 
matching yields the following dispersion relation for a propagation constant of a mode in waveguide with anisotropic 
core: 

[Jl(K TE R) + Ki(n cl R)] [e {l Jl(K TM R) + e cl Kl(n cl R)] = + (-^- + (7) 

U C \K TE K cl J \K TM K cl J 

where L^kR) — L' m (nR) /[nL m (KR)]. The two terms in the left-hand side of the equation correspond to the 
contributions from TE and TM modes respectively. As it is follows from Eq.© the "pure" TM and TE modes are 
only possible when (i) m = 0, or (ii) e c \ — > — oo. The latter case corresponds to perfectly conducting metallic walls 
described in Section 1. Solutions of Eq.Q can be separated into two fundamentally different groups: the ones with 
k (te\tm) ^ describe volume modes, while the ones with k 2 te i tm -> < correspond to surface waves. 

It is possible to dramatically simplify Eq.© for the case of waveguides with metallic walls. At optical or infrared 
frequencies, the permittivity of metals is dominated by plasma-like response of their free electrons^. As a result, this 
permittivity is negative and |e c i| 3> 1. A straightforward Taylor expansion of Eq.JTI) yields 



,TE 



ife< > 



ck { z 0)2 R yf^ci 
c4 ) *J m ($ ) E R) ( fcgm^j \ 1 
. uk^RJ^lR) \ e^n^R 2 ) ^ 



where the superscript (0) denotes the mode parameters in a waveguide with perfectly conducting walls. Note that 
similar to planar waveguides^, finite value of permittivity of the waveguide wall leads to a mode expansion into the 
cladding region. 
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FIG. 4: TMii(HEu) volume modes in waveguides with anisotropic cores. Permittivities of core material correspond to Si-Ag 
composite described in the text. Ao = 1.2/im. Panel (a) corresponds to perfectly conducting waveguide walls, panel (b) shows 
the mode in a waveguide with Ag cladding. Lines represent exact results [Eq. Q]; dots correspond to perturbative Eq.JSJ 




FIG. 5: TM02 mode in anisotropic fiber with free-standing anisotropic core (e c i = 1) (a); SPP mode propagating at Ag-Air 
interface (b); Ao = 1.2/im. Note the fundamental structural difference between volume and surface modes. Also note that 
volume mode has better confinement than its surface counterpart. 



Besides affecting the mode propagation constant, cladding material in fiber systems also affects the mode structure. 
In fact, the m > 1 mode of a cylindrical waveguide with real (metal or dielectric) walls can be represented as a linear 
combination of TE and TM waves, known as HE or EH waves^. In particular, the HE wave can be represented as a 
combination of the TM mode with an admix of the TE mode: E z oc Jmii^TMf) exp(im<p), H z cx aJm^KTEf) exp(im0), 
with the admix parameter a given by: 

iwfe z me|i(e_L - e ci ) J m (n TM R) 

ce ± K 2 cl K*, M R jUk te R) + KL( Kcl R) 

Note that the effect of mode structure modification is unique to fiber-geometries and is not observed in planar 
waveguides. Our calculations show that for Ag walls, this admix in the HE\\ mode is below 2%. 

As noted before, the TE (EH) modes are (almost) unaffected by the material anisotropy. Therefore the properties 
of these waves are identical to the properties of TE modes in waveguides with isotropic core. 

In Fig0]we illustrate a propagating volume mode in the Ag-Si system described above for Ao = 1.2/xm. We provide 
a comparison of the mode structure in waveguides with perfectly conducting and Ag walls. It is clearly seen that for 
the silver waveguide, the mode structure is well-described by the perturbative result [Eq.®]. 

Finally, use the exact dispersion relation [Eq.Q)] to calculate the modes in the metamaterial fiber without cladding, 
and compare these modes to the surface polariton mode on the Ag nanowiraSS. Results of these simulations are shown 
in FigJSJ It is clearly seen that the structure of the surface mode (localized at core-cladding boundary) is fundamentally 
different from the one in the volume modes. It is also seen that surface waves have weaker confinement than their 
volume counterparts. 
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V. CONCLUSIONS 



In conclusion, we have presented an analytical description of light propagation in waveguides with strongly 
anisotropic dielectric cores. We demonstrated that these systems support propagating modes even when the waveg- 
uide radius is as small as 1/15 of free-space wavelength. We have analyzed the effect of material microstructure 
and waveguide cladding on mode propagation and structure, and suggested a practical realization of subwavelength 
structure for telecom frequencies. 

Finally, we note that our analysis can be easily generalized for the case of non-circular waveguides, for different 
classes of dielectric, plasmonic, or polar materials, and for different frequency ranges, from optics to IR to THz 18 . 

This research was partially supported by the General Research Fund, Oregon State University. 
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